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I. Technical Summar 


l. Background and Accomplishments 

The research performed under this grant had as its objective the 
development of a ray-optical theory for the behavior of eigenvalues and 
eigenmodes in unstable laser resonators. The detailed program was de¬ 
termined za co isultation with Drs. J. Hanlon, R. Butts, and G. Dente at 
the Air For^c Weapons Laboratory, Albuquerque, New Mexico, who served 
successively as technical project monitors. The investigation proceeded 
along two principal lines. 

1. Waveguide Analysis 

Further exploration of the waveguide analysis of Chen and Felsen* 

2 

and Santana and Felsen was undertaken to explain, on a coupled waveguide 
mode basis, the structure of the eigenmode charts obtained by the con¬ 
ventional "brute-force 11 numerical solution of the resonator integral equation. 
Here, the resonator is regarded as a waveguide perpendicular to its axis, 
with the mirrors representing the waveguide boundaries (Fig. 1). Coupling 
between selected modes in this highly overmoded structure is introduced 
by the mirror edges which represent waveguide discontinuities. The scatter¬ 
ing from these discontinuities is taken into account by a ray-optical analysis 

3 

incorporating the constructs of the geometrical theory of diffraction (GTD). 

These edge-diffracted fields are then converted into modal coupling coef¬ 
ficients to permit the formulation of a waveguide mode resonance equation 
("transverse resonance" equation) that specifies the eigenvalues of the composite 
resonator modes. Thus, the method involves a combination of waveguide 
mode and ray-optical techniques. 

The method had previously been applied successfully to sharp-edged, 

2 

aligned, symmetrical strip and circular mirror configurations. Its ex- 

4 

tension to misaligned sharp-edged strip mirrors (Fig. 2) and to aligned strip 
or circular mirrors with rounded edges (Fig. 3) was undertaken under this 
grant.^ The results for misalignment were compared to numerical data 
based on the resonator integral equation, and were found to be interpretable 
in terms of coupling between waveguide modes with even and odd symmetry 
with respect to the resonator axis. The coupling disappears when the mirrors 
are perfectly aligned. New results for very small misalignment, inadequately 
explored in the literature, were found from a careful study of the transverse 









resonance equation for different ranges of equivalent Fresnel numbers, 

The results obtained for round-edged mirrors were entirely new.^ 
Of particular interest was the demonstrated existence of an azimuthally 
symmetric low-loss detached mode for spherical mirror structures; mode 
detachment here is not observed in the sharp-edged mirror case. For 
strip mirrors (two-dimensional case), edge rounding was shown to intro¬ 
duce detachment of the azimuthally symmetric mode at lower values of N 
and increased mode separation (see Appendix A). The waveguide method 
was shown to be effective for simple calculation of the eigenvalues, and 


for cogent interpreation of the complicated eigenmode charts. This is due 
to the fact that coupling between only two systematically selected modes is 
adequate to explain the resonances. However, the more sensitive eigen¬ 
mode fields on the mirrors are not well enough expressed in terms of this 
simple model; coupling to many modes would be required to yield adequate 
mode shapes. For this reason, the waveguide method is not regarded as 
an attractive alternative for computation of unstable resonator modal fields. 


-Optical Analvsis 


By an alternative approach, an eigenmode in the unstable resonator 
may be regarded as being established by self-consistent treatment of the 
mirror-edge-diffracted ray-optical fields that find their way back into the 
resonator by direct and multiply reflected paths between the mirrors. This 
direct ray tracing method, performed according to the rules of GTD, is 
hampered by the existence of ray fields with very many reflections, there¬ 
by making direct summation over these fields impractical and inaccurate, 

A breakthrough in our analysis of this problem has been the ability to deal 
collectively with the higher order reflected rays. ' Thus, the resonant field 
is represented in terms of a selected number of lower-order reflected rays 
plus a collective ray that accounts for the cumulative effect of all higher 
order reflections. This procedure was found to yielc not only a numerically 
attractive and accurate alternative to the "brute-force 11 method of integrating 
the resonator integral equation, but it also furnishes some penetrating in¬ 
sights into the behavior of the unstable resonator in purely ray-optical 
terms, without recourse to the resonator integral equation. 


This study was motivated by the asymptotic theory for the resonator 

8 9 10 

integral equation developed by Horwitz and applied by him 7 and others 
to various resonator configurations. Here, it was shown that the solution 
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of the eigenvalue problem could be reduced to solution of a polynomial 
equation, and that the corresponding eigenmode fields can be represented 
in terms of a set of certain wave functions. The method, which was found 
to yield remarkable agreement with eigenvalues and eigenmode shapes ob¬ 
tained by the direct numerical procedure, is based on a postulated form of 
the solution (ansatz). While the postulates were qualitatively identified 
with various ray-optical phenomena, no quantitative correspondence and 
interpretation was established. 

Our ray-optical method provides a complete and quantitative ray- 
optical theory of the unstable resonator. The results of the Horwitz method 
emerge from our analysis as a special case, and each element in the 
Horwitz ansatz is derived deductively. To test our method, it was first 
applied (with complete success) to symmetric, aligned, sharp-edged strip 1 
and spherical mirror** structures. It was also applied to the round-edged 
mirror case. *^ From the understanding of the physical phenomenology 
gained thereby, it appeared feasible to attack more systematically other 
resonator configurations that depart from the idealized shapes investigated 
originally. A first step in this direction was an application of the ray method 
to half-symmetric resonators with internal axicon (HSURIA), *^ which pro¬ 
vided an understanding of the effect of the internal axicon reflector, and of 
the influence of shielding of the axicon tip. 

B. Summary of the Ray Method 

As noted in Section A2, a key feature in our procedure has been the 

ability to show that ray-optical phenomena descriptive of high-order multiple 

reflections can be treated collectively . Thus, our approach involves the 

simultaneous use of individually tracked ray fields as well as collective ray 

fields that resemble modal congruences. Modal ray congruences describe 

2 

self-consistent fields in resonators with infinite mirrors, and they are 
exemplified in a special case by the two self-replicating cylindrical wave 
fronts of Siegmann. More generally, these wave fronts, and the associated 
normal ray congruences, are generated by elliptical caustics when the 
infinite mirrors ha ,r e hyperbolic shape. 


A second key feature has been our ability to construct the resonance 
equation for the eigenvalues and eigenmodes by invoking self-consistency 
of the modal fields without recourse to the integral equation. Thus, both 











the fields per se and the resonance condition have been formulated entirely 
within a deductive ray-optical framework. 

Our procedure for the symmetrical strip resonator involves the 
following steps 

1. Determining the collective (modal) ray fields (self-replicating 
wave fronts and ray congruences) that describe the fields in a resonator 
with infinite hyperbolic mirrors. This incident modal field, which can be 
found by direct ray-optical techniques or by the WKB method, provides the 
excitation. Its precise form is determined by the elliptic caustic which is 
characterized by the initially assumed caustic parameter, p (Figs. 4(a) 
and (b)). 

2. Accounting for edge diffraction, when the mirrors are finite, 
by postulating equivalent edge line currents with non-isotropic radiation 
pattern (Fig. 5), The pattern function f( 0 , 6.) in the direction 0 is deter¬ 
mined from the known solution of edge diffraction by a half-plane when an 
incident plane wave arrives at the angle 6.. The local modeling of the 
curved mirror near the edge by a plane surface and of the modal ray field 
incident on the edge by a plane wave is legitimized by the rules of GTD. 

The starting amplitude of the field on each edge-diffracted ray thus depends 
on 0. which, in turn, depends on the caustic parameter p. 

3. Tracking of edge-diffracted ray fields into the resonator. This 
tracking takes place via the well-known rules of geometrical optics, relying 
on optical path lengths to establish phase (Fig. 6(a)) and on conservation of 
energy in a ray tube to establish amplitude (Fig. 6(b)). Various species of 
multiply reflected rays are identified in Fig. 6(a). Some of these are re¬ 
flected out before they reach the resonator axis while others pass through 
the resonator from one side to the other. 

4. Summing of edge-diffracted ray fields with many reflections in¬ 
to a collective closed form. These collective rays appear to emanate from 

the foci of the hyperbolic mirror system and thus resemble the self-replicating 
modal ray fields that provided the initial excitation (Fig. 7), 

5. Determining the total singly diffracted ray-optical field. The 
field due to single diffraction at each edge, as observed at a point within 
the resonator, consists of 2N individually reflected ray fields plus the field 
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carried along the collective ray, which accounts for all ray fields having 
undergone more than N reflections (Fig, 8). The choice of N is not critical 
but must be large enough to legitimize the collective treatment. 

6. Determining the ray-optical field due to double diffraction. When 
the observation point is chosen at an edge, the (2N+1) singly diffracted ray 
fields in Fig. 8 provide the excitation for double diffraction. The resulting 
doubly diffracted ray fields are tracked through the resonator in precisely 
the same fashion as the singly diffracted fields. However, the excitation 
amplitude of each ray species now depends on the angle of incidence 8^ of 
each of the (2N+1) singly diffracted rays, thereby requiring the pattern 
function f( 8 , 8^) instead of f( 8, 8^) (Fig. 9). 

7. Establishing the resonance condition. Resonant solutions for 
the unstable open resonator require self-sustaining fields in the absence of 
excitation. Thus, the field incident on an edge after a certain number of 
diffractions must be identical with the field incident on the same edge before 
the last diffraction took place. Therefore, with reference to the above, the 
singly diffracted fields incident on an edge must be the same as the doubly 
diffracted fields (Figs. 8 and 9). Imposition of this condition yields the 
equation for the resonant eigenvalues. These eigenvalues are directly re¬ 
lated to the caustic parameter p in Fig. 4a. In essence, the resonant p 
values adjust the initially assumed incident ray congruences so that the re¬ 
sulting field in the finite mirror structure is self-resonant. 

8. Determining the eigenmode fields. These are given by the (ZN+1) 
ray-optical fields when the resonant eigenvalues are inserted for the caustic 
parameter p that appears in each of the ray field expressions. 

The details of the above-described formulation have been given in 
references 7 and 13. It has also been shown there how each of the wave 

g 

functions and recursion relations in the Horwitz ansatz can be interpreted 
ray-optically, and how the ray-optically derived resonance condition reduces 
to that of Horwitz after certain approximations. The same correspondence 
between the ray-optical approach and the procedure of Butts and Avizonis^ 
for symmetric circular mirror configurations has also been demonstrated, * 
Here, the simple ray-optical method requires correction near the resonator 
axis where edge diffracted fields due to the circular rim accumulate. 
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The successful ray-optical solution of the ideal symmetric strip and 
circular mirror resonators has demonstrated two features with major im¬ 
plications for more general problems: 

1. The ability to construct the field within the resonator by local 
tracking along geometric-optical ray paths, with the possibility of summing 
ray fields with many reflections into a collective form. 

2. The ability to construct a resonance equation by invoking self- 
consistency on the ray fields before and after a successive diffraction event. 

Feature 1 implies that modifications of the ideal geometry can be 
accounted for by corresponding modifications of the appropriate ray fields. 

This is particularly relevant when these modifications occur locally (as 
when the medium is inhomogeneous, or when there are local perturbations 
on the mirrors) since the rays describe locally tracked fields. The ray 
tracking also clarifies obscuration and blocking effects introduced by axicons 
in annular resonator configurations, by scraper mirrors or by coupling 
apertures. 

Feature 2 implies that once the ray structure and the ray fields have 
been determined, these fields can be used directly to construct the resonance 
equation without the need for going through the resonator integral equation. 

In fact, by appropriate categorization of various ray species, one may ascertain 
the influence of each on the resonance behavior. 


These aspects are illustrated for the HSURIA (half symmetric unstable 
resonator with internal axicon). Various possible ray species are depicted 
in Fig. 10. They are ordered like those in the equivalent symmetric circu¬ 
lar mirror configuration, which can be constructed from the HSURIA by 
imaging and unfolding. Due to the blocking effect of the axicon, some of 
the rays of the unfolded structure are eliminated. The categorization into 
existing and blocked rays can be undertaken directly in the symmetric 
resonator if the position of various rays relative to the axicon tip location 
is observed (Figs. 11 and 12). The resulting resonance condition is the 
same as for the symmetric unfolded resonator provided that one eliminates 
therefrom the blocked rays. However, fields on rays reflected near the 
axicon tip should be modified to account for the blocking ( ,f reflected-ray 








boundary" effect). Tip diffracted rays could be included to refine the 
resonance condition. 

D. Conclusion 

The ray method has been developed sufficiently to warrant applica¬ 
tion to various problem areas related directly to high-power laser systems 
now under development, and to provide alternatives to, and independent 
checks of, the purely numerical procedures now in use. Some topics for 
investigation have been suggested by us, which address remaining questions 
in the conventional symmetrical mirror unstable resonator , and also in annular 
resonators with conical elements, which are gaining increased importance. 
However, these studies, which include effects of scraper mirrors, coupling 
apertures in injection-locked resonators, medium mhomogeneities, etc., 
could not be performed within the grant period. 
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Resonator 

axis 



parallel plane 
regions 


Fig. 1. Waveguide model of the unstable resonator with hyperbolic 

mirrors. A waveguide mode in the infinite mirror structure 
(solid plus dashed curves) is described by an elliptic caustic 
with foci at i d and caustic parameter p. Truncation of the 
mirrors (solid curves only) causes waveguide mode coupling 
and reflection. To determine the mode coupling and reflection 
coefficients, the local regions near the edges may be regarded 
either as equivalent parallel plane regions (for small Fresnel 
numbers), or as individual equivalent half-plane scattering 
centers (for large Fresnel numbers). Resonance defining the 
resonator eigenmodes is established by invoking a "transverse 
resonance" condition on the waveguide modes, which ensures 
self-consistent reflection of modal fields traveling to the 
right (f?) and left (R). 
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a) misaligned resonator b) equivalent resonator with unequal mir- 

(no symmetry with ror length (symmetry with respect to 

respect to axis) axis); n] and n 2 are hyperbolic sur¬ 

faces in an elliptic coordinate sys¬ 
tem. (3 g and (3 are the angles 
made by the modal ray at the edges 
B and C, respectively. 

Fig. Z - Resonator with misaligned circular strip mirrors 











Edge termination 



Fig. 3. Resonator with rounded edges. At y = y, a section of 

waveguide with circular cylindrical boundaries of radius 
r£ is joined smoothly to a cylindrical resonator which 
otherwise would be sharp-edged. 
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Fig. 4(a). Ray congruences for self-consistent field between infinite, 

hyperbolic strip mirrors: caustic for 2kL^nrr. The caustic 
parameter p is related to the ellipse u= « 
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Fig. 7. Singly diffracted ray fields G^“(x) reaching an observation point x 

on the lower mirror after (n-l) reflections; the incident rays (shown 
dasned) oelong to the modal field between infinite hyperbolic mirrors. 
G-(x) represents the collective ray field for n> N. 
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right-hand edge due to primary diffraction 


Total ray field at lower 
of incident field. 
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ijectories to A for various ray species, dealing directly with the 
oided structure. The dashed ray n = -1 and others :n this categ or\ 
t shown) 3 re eliminated :v the axitor.. 



















axicon tip 



Fig. 11(b). Equivalent construction of the domain in Fig. 11(a) (shown 

shaded) directly in the unfolded configuration, noting placement 
of the axicon tip. 
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RAY-OPTICAL CALCULATION OF EIGENMODE BEHAVIOR OF UNSTABLE 


USER RESONATORS WITH ROUNDED EDGES 

C. Santana 

Institute de Pesquisas Espaciais - INPE 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico - CNPq 
Sao Jose dos Campos, SP, Brazil. 


ABSTRACT 


The eigenmode resonance equation previously developed by the ray-optical 
method of analysis for unstable, symmetric, bare resonators with strip and 
spherical mirrors is employed here to establish the basic properties of resonators 
with rounded mirror edges. The numerical results obtained in this manner are far 
more comprehensive than those presented earlier from a waveguide model approach 
to the resonator problem. As in this earlier study, edge rounding is found to 
improve mode separation and to shift the onset of low-loss mode detachment to 
lower Fresnel numbers. It is also found that a very simple closed form expression 
predicts quite accurately the low-loss mode behavior after it has detached from 
the remaining mode spectrum. 


- 30 - 









I. Introduction and Summary 

Edge rounding of strip mirrors in unstable symmetric laser resonators 
initiates detachment of the low-loss eigenmode at small values of the equivalent 
Fresnel number. For the spherical mirror geometry, rounded edges produce a 
low-loss detached mode that is not found in the sharp-edged configuration. Mode 
discrimination is improved in both cases. Because of these attractive features 
which have been pointed out in a previous publication 1 , it is desirable to 
generate detailed eigenmode charts that exhibit the dependence of the resonance 
on the edge radius. The resonator edges are modeled by an abrupt change in the 
mirror curvature as shown in Fig. 1. The termination of the added mirror section 
is assumed to be constructed so as to render its effect negligible. 

Some preliminary results obtained via the waveguide approach to the resonator 
problem have already been presented 1 . The waveguide resonance equation is 
convenient for predicting the physical mechanism responsible for the resonance 
and for dealing with the lovest-order detached mode, but it is less easily 
implemented for determination of the behavior of the higher order modes. An 
alternative view of the resonant wave process is obtained by a ray-optical 
analysis 11 " 4 which leads to a resonance equation that is equally applicable to 
the lowest order and higher-order modes. The ray-optical resonance equation for 
round-edged mirrors 4 - is employed for the generation of the data in this 
communication. 

II. Strip Resonators 
A. Resonance Equation 

From the analysis in Ref. 2, the eigenvalues X for the even and odd symmetric 
modes in unstable strip resonators are given by the polynomial equation 
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The power loss per transit is related to the eigenvalue \ and the magnification 


PL * 100 (1 - U| 2 /M) . 



For the N values for which the resonator loss approaches the geometrical 
optics value (corresponding to j\| - 1), an approximate equation for the eigenvalue 
of the even symmetric mode, which is the lowest-loss mode 5 , can be found by 
iteration as in Ref. 3: 


X = 1 


2e2 exp(i2*rN T eq ) 


l-2e2^exp (i^N ) 


provided that 

£n(47rN e q) /£nM << N << S,nM/2 | e 2 ! . (6) 

This simplified equation developed from the rav-optically derived polynomial 
resonance equation is in agr ~ent with that obtained from the waveguide analysis. 
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B. Numerical Results 

The full resonance equation (1) has been solved numerically for the 

eigenvalues of an unstable strip resonator with rounded edges having a 

magnification factor M * 2.9. Figure 2 shows the magnitude of \ versus N for 

eq 

an edge rounding factor r/ri = 100. It is seen that, in conformity with the 
results of Ref. 1, the essential effect of edge rounding in strip resonators is 
the formation of a low-loss detached mode at lower Fresnel numbers than for the 
sharp edge case. In a sharp edged strip resonator with M * 2.9, mode detachment 
accurs only for N e ^ > 9, 5 in contrast to > 4 as seen from Fig. 2. A 
comparison of Fig. 2 with the plot of j\j for a sharp edged resonator with the 
same parameter values^ reveals that edge rounding also enhances mode 
discrimination. 

Figure 3 shows a plot of j\| versus N ecj for the same resonator as in Fig. 2, 
but with r/ri * 30. We note that as the edge effect is deemphasized by decreasing 
r/r 2 , the geometrical optics value is attained for smaller N e q values and the 
separation between the low-loss detached mode and the next higher order modes is 
increased. 

The simplified eq. (5) can be utilized for the calculation of X in the N e ^ 
region where jxj 2 : 1, which generally corresponds to the detached mode region. In 
fact, for N e q > 2.5, the j X| values for the lowest-loss modes of Figs. 2 and 3 
if calculated via eq. (5), deviates negligibly from the values calculated by the 
full resonance equation. Equation (5) is also in excellent agreement with the 
power loss calculations effected by modal analysis 1 . 

In the numerical comput-* ions, it has been found that the parameter N, 
which characterizes the number of non-modal rav fields can be set in accordance 
with Horwitz’s criterion * 250 N e q. 
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For plotting the eigenvalue cui.es, the roots of the full resonance equation 
(1) were found at each N e q value, in steps of 0.1, except in the vicinity of 
those N e q values where mode detachment was to be ascertained. There, the step size 
was reduced to 0.01 and the detailed evolution of the complex values of X was 
determined in order to show which modes evolved into which. It is found that if 
the step size is not reduced, one can easily arrive at a false conclusion about 
the onset of mode detachment. 

III. Spherical Mirrors Resonators 
A. Resonance Equation 


The polynomial equation for the eigenvalues * of an unstable resonator with 
round-edged circular hyperbolic mirrors can be found from the results of Ref. 2: 


aM m - m N [ i + - (-i)™ 


L 


t=i 


(1 - M ~ K ) 2 F, (1) 

t ,m 


m. 


2irN i m 

(-1 ) m exp(i2rN eq > | j 


* 0 


(7) 


where 


j m i(r/r 2 - 1) M - 1 


4rN 


eq 


F (x) * 2 exp , 
t, m 


M * 1 
i2 7 rN 




l-M 




. j 


j J m <V [ 1 * 3 * 2 a " : ' ' - 1 J ..I <V * M '' L 


3 + x- M 


•21 1 


(1 - x 2 sf‘V 


( 8 ) 


(9) 


and the remaining terms are in conformity with the definitions of Ref. 4, 

The power loss per trar is related to ;•, and M via 

PL = 100 [ 1 - (!X : /M) 2 ] . (10) 
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As in the cylindrical mirror case an approximate equation for the 
eigenvalues, valid when the 1 ss approaches the geometrical optics limit \=l 
for the azimuthally symmetric modes (m=0), can be found by iteration of eq. (7): 


6exp(i2TrN ) 

X = 1-——-S3- 

1 - 0 N exp ( i2 ttN ) 
eq 

with the restriction 

£n(4rN )/£nM N << fcnM/6 

eq 


( 11 ) 


( 12 ) 


B. Numerical Results 

To assess the effects of edge rounding in spherical mirror unstable 

resonators, the full resonance equation (7) was solved numerically for various 

resonator parameters. Figures 4 and 5 show plots of ) X | versus N for a resonator 

with M®2 and r/r£ - 100* The most obvious difference between these plots and their 

correspondents for a sharp edged resonator 5 is the appearance of a detached mode 

for N > 59, thereby confirming the prediction made by modal analysis, that 

mode detachment does take place in three-dimensional resonators with rounded 

edges 1 . One may also observe that the eigenvalue approaches the geometrical optics 

limit I\| = 1 faster than in the sharp-edged configuration. Another characteristic 

of edge rounding is illustrated in Fig. 6, which is a plot of !\ versus N for 

the same resonator parameters as in Figs. 4 and 5, but this time for the lowest 

order asymmetric (m a l) modes. Figure 6 should be compared with Fig. 5 of Ref. 6 

to conclude that edge rounding improves the discrimination between the dominant 

modes for m*0 and m*l, even in the N range where mode detachment for m=Q has 

eq 

not yet been attained. The ' ^gner losses achieved by edge rounding for the m*l 
modes are welcome since these modes are influential in beam deterioration for not 
having a maximum at the center. All of these features caused by edge rounding are 
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in agreement with the predictions made by modal analysis, 1 where it was 

anticipated that the eigenvalue for the m*l modes should approach the geometrical 

optics limit •>, = 1/M onlv for very large N values. 

eq 

Depending on the edge rounding parameter r/r^, mode detachment for three- 

-dimensional resonators may occur for low Fresnel numbers, as is the case 

illustrated in Fig. 7, where |\| is plotted versus N for a resonator with 

M = 2.5 and r/r? =* 16.7. Here, detachment starts for N 6. 

z ’ eq 

The simplified equation (10) gives good accuracy in the detached mode region 
and is also in very good agreement with the power loss per transit calculations 
by modal analysis. 4 
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FIGURE CAPTIONS 


Physical configuration. For strip or circular mirror resonators, the 

rounded edge is modeled by a circular section with radius r 2 , while 

the radius of curvature of the main mirror is r. 

Modulus of \ versus N for a strip mirror resonator with M * 2.9 and 

eq 

r/r 2 = 100. Even symmetric modes. 

Modulus of a versus N for a strip mirror resonator with M * 2.9 and 

eq 

/o 35 30. Even symmetric modes. 

Modulus of a versus N for a circular mirror resonator with M * 2 an 

eq 

r/r 2 * 100. Azimuthallv symmetric modes. 

Modulus of a versus N for a circular mirror resonator with M * 2 an 

eq 

r/r 2 * 100. Azimuthally symmetric modes (m*0). The stationary roots 

have been suppressed. 

Modulus of a versus N for a circular mirror resonator with M * 2 an 

eq 


r/r 2 ■ 100. First asymmetric modes (m=l). 


Modulus of * versus N for a circular mirror resonator with M ■ 2.5 

eq 

and r/r 2 * 18.0. Azimuthally symmetric modes (m*0). 
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